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INTRODUCTION 

Tne  optimum  design  of  aircraft  structures  and  the  maximum  exploitation 
o.  the  auvantageous  properties  of  advanced  fiber  reinforced  materials  re- 
quire  detailed  knowledge  of  the  stress  fields  within  structural  elements. 

Of  particular  in?>ortance  are  stress  concentrations,  for  it  has  been  indicated 
that  they  are  of  a  different  nature  in  composites  than  in  conventional 
Mterials.  Under  reference  (a),  the  NAVAIRDEVCEN  (Naval  Air  Development 
Center)  has  conducted  a  program  to  evaluate  the  stress  concentrations  in 
composite  structures, 

Tne  most  basic  form  of  stress  concentration  in  a  composite  laminate  is 
-he  occurrence  of  interlaminar  shear  and  normal  stresses  near  free  edges, 
accompanied  by  variations  in  in-plane  stresses.  This  stress  concentration, 

effect,  is  caused  by  the  heterogeneity  of  the 
Inmate,  ai^  is,  therefore,  a  characteristic  more  of  the  material  than  of 
the  geometrical  configuration  of  a  structural  element.  The  resulting  per- 

within  a  distance  frota  the  edge  on  the 
order  of  the  laminate  thickness,  but  may  be  sufficient  to  cause  delamination 
and  subsequent  failure. 

Since  conventional  laminate  theory  does  not  account  for  the  edge  effect, 
an  alternate  method  is  needed  to  analyze  this  localized  phenomenon.  Several 
approaches  to  the  special  case  of  a  straight  edge,  using  a  variety  of 
analytical  and  numerical  methods,  have  been  published.  (References  (b) 
through  (f)).  However,  the  methods  used  are  either  two-dimensional  or  are 
applicable  only  to  straight  edges.  The  effort  at  the  NAVAIRDEVCEN  has  been 
the  development,  using  NASTRAN,  of  a  three-dimensional  finite  element 

and  analysis  technique  which  can  be  applied  to  edges  of  arbitrary 
geometiy  including  holes  and  notches,  and  to  bonded  joints.  The  technique 
developed,  and  its  application  to  the  problem  of  a  four-ply  ±45“  laminate 
in  tension  are  described  in  this  report.  The  results  of  the  analysis  of 
tnis  elementary  case  are  essentially  identical  to  those  obtained  by  other 
^vestigators  and  indicate  that  the  interlaminar  shear  stress  between  the 
+45  and  -45  plies  may  be  sufficient  to  cause  delamination. 

PRELIMINARY  ANALYSIS 

Because  of  the  availability  of  solutions  by  other  methods,  the  finite 
element  procedure  has  been  first  applied  to  the  free  edge  of  a  four-ply,  ±45" 

tension,  as  shown  in  Figure  1.  It  is  postulLL  that 
xn  the  middle  region  of  the  laminate,  away  from  the  areas  of  load  introduction. 

transverse  displacement 

in  the  longitudinal  ( ?C.  )  direction  (parallel  to  the  edge);  only  the  varia- 

addL°e  displacement  m  )  is  non-zero.  In 

addition,  by  introducing  compatibility  relations  from  reference  (g)  and 
setting  all  partial  derivatives  with  respect  to  X  equal  to  zero,  it  is 
concluded  that  is  a  linear  function  of  y  and  z  of  the  forra 

Ay  +  B2  +  C.  However,  since  the  laminate  and  loading  are  symmetric,  the 
linear  ternas  must  be  zero,  and  therefore,  is  constant  over  the  cross 

section.  Aiso  oecause  of  symmetry,  there  are  no  shear  stresses  on  the 
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mid^piane  of  the  laminate j,  nor  is  there  any  dispiacement  nomal  to  it* 

The  most  important  consequence  of  the  lack  of  variation  in  the  longitudinal 
direction  and  of  the  problem  symmetry  is  that  only  two  piles  of  a  thin 
slice  normal  to  the  edge  need  to  be  modelled  and  analyzed^  In  addition^ 
since  the  edge  effect  decays  by  97  percent  within  a  distance  equivalent  to 
the  laminate  thickness  (four  plies)  ,  the  region  to  be  analyzed  extends  only 
six  ply  thicknesses  from  the  edge,  to  where  the  stress  state  can  be  deter¬ 
mined  by  conventional  laminate  theory The  region  to  be  analyzed,  and  the 
applicable  stress ^  strain  and  displacement  conditiotis  are  shown  in  Figure  2. 

FINITE  ELEMENT  MODELLING 

The  grid  point  array  of  the  firdte  element  model  is  shown  in  Figure  3, 

To  facilitate  modelling  the  orthotropic  material,  the  mesh  is  square  in 
planform,  subdivided  to  one-eighth  ply  thickness  near  the  edge,  where  stress 
gradients  are  greatest^  A  single  ply  of  the  laminate  is  represented  by 
four  planes  of  grid  points,  as  indicated  by  the  shaded  region  in  Figure  3. 

The  grid  points  are  interconnected  by  one-  and  two-dimensional  elements  -- 
rods,  shear  panels  and  quadrilateral  membranes  --  which  together  siitiulate 
the  three-dimensional  elastic  behavior  of  the  orthotropic  materials 

An  enlarged  view  of  the  shaded  region  of  Figure  3  is  shown  in  Figure  4, 
The  boron-epoxy  material  is  modelled  as  two  layers  of  epoxy  separated  by  a 
fibrous  layer  whose  thickness  is  equal  to  the  fiber  diemeter.  The  properties 
of  the  fibrous  layer  (Table  I)  are  derived  from  those  of  the  unidirectional 
material  and  the  epoxy  by  means  of  a  three-dimensional  laminate  analysis. 
Since  each  of  the  layers  is  orthrotropic,  its  elastic  behavior  can  be  de¬ 
scribed  by  nine  independent  stiffness  coefficients  (Figure  5).  Ttie  stiff¬ 
ness  represented  by  each  coefficient  is  provided  in  the  finite  element  model 
by  one  or  more  of  the  elements.  For  example,  the  interlaminar  membranes, 
which  are  aligned  pai-aliel  and  normal  to  the  fiber  directions,  contribute  to 
the  stiffness  represented  by  the  coefficients  €5^3  and  C23.  The  transverse 
normal  stiffness,  C33,  and  the  Interlaminar  shear  stiffnesses,  and  C55, 
are  provided  by  the  rods  and  shear  panels. 

LOADING  AND  CONSTRAINTS 

Tile  knov;n  stress  and  displacement  conditions  applicable  to  the  analysis^ 
have  been  shown  in  Figure  2  and  the  constraints  imposed  on  the  model  to 
simulate  the  strain/displacement  conditions  are  listed  in  Table  II.  Zero 
dispiacement  across  the  laminate  mid-pIanc,  the  lack  of  variation  of  dis¬ 
placement  in  the  longitudinal  direction,  and  the  uniform  norraal  strain, 

^  ,  are  enforced  by  vaeans  and  multi-point  constraints.  At 

the  interior  end  of  the  model  ^  K  ^  )>  additional  constraints 

{$<0/8^-  6U/6Z-  6V/5Z^  0)  are  applied  to  simulate  the  laminate 

theory  soIutiori>  No  forces  are  input.  Instead,  the  unknown  grid  point 
forces  are  generated  by  NASTRAN  as  constraint  forces  in  response  to  a 
sufficient  irumber  of  known  displacement  conditions. 
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STRESS  RECOVERY 

To  calculate  the  stresses  in  the  individual  fibrous  and  epoxy  layers, 
the  elaaent  stresses  are  transformed  as  necessary  into  a  common  coordinate 
systen  (laminate  coordinates)  and  added  together.  For  example,  if  the 
longitudinal  normal  stress ,  (J^  is  expressed  as 

^  *  C,(  €x  +  C|£  Cy  C|5  +  Ctg 

the  first  three  terms  are  equivalent  to  the  average  of  the  stresses.  Ox  , 
from  the  upper  and  lower  in-plane  membranes  and  the  third  term  is  the  sum 
of  the  stresses  (transformed)  of  the  interlaminar  membranes.  The  stresses 
obtained  by  this  method  are  applicable  for  the  middle  of  each  layer.  In¬ 
plane  ply  stresses  are  obtained  as  the  sxaa  of  the  fibrous  and  epoxy  layer 
stresses,  each  weighted  in  proportion  to  its  thickness. 


RESULTS  AND  DISCUSSION 

The  ^-plane  ply  stresses,  Ox  and  ,  and  the  interlaminar  shear 

stress,  the  adjoining  epoxy  layers  of  the  two  plies  are  plotted 

versus  distance  from  the  laminate  edge  in  Figure  6  and  are  essentially  iden¬ 
tical  to  those  obtained  by  other  investigators.  The  remaining  stresses, 

0/  »  and'*)Cy^  ,  which  arise  primarily  from  material  dissimilarities 

betx,reen  the  fibrous  and  epoxy  layers,  are  two  orders  of  magnitude  smaller. 

At  the  laminate  edge,  where  the  in-plane  shear  stress , ,  is  zero,  the 
longitudinal  tensile  stress,  CTx  »  is  the  same  as  would  occur  in  unidirec¬ 
tional  material  loaded  to  the  same  normal  strain,  ,  at  45  degrees  to 

its  fiber  direction,  and  is  20  percent  lower  than  in  the  interior.  In  order 
to  maintain  equilibrium,  the  rise  of  the  in-plane  shear, ^xy»  its 

asyt^totic  approach  to  the  value  predicted  by  laminate  theory  are  accompanied 
by  a  decay  of  the  interlaminar  shear  ,  from  its  maximum  value  at  the 

edge.  For  each  of  these  stresses,  the  influence  of  the  free  edge  decays  by 
97  percent  within  four  ply  thicknesses.  From  an  attenqjt  to  fit  the  data 
obtained  from  the  finite- element  analysis  to  analytical  forms,  it  has  been 
found  that  the  variation  of  the  major  stresses  in  the  direction  normal  to 
the  edge  can  be  very  precisely  described  by  simple  exponential  functions, 
as  shown  below; 

a;  -  -0^)0 
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where  h^is  the ^ ply  thickness,  and  the  superscripted  bar  and  zero  denote 
vuiues  in  the  interior  and  at  the  edge,  respectively o 


Ti-xe  magnitude  of  the  interlaminar  shear  stress  at  the  edge,  about  equal 
to  the  inplane  shear  predicted  by  laminate  theory  for  the  interior,  is  most 
important,  as  the  occurrence  of  so  high  a  stress  in  the  matrix  may  cause 
delamination  and  redistribution  of  stress  leading  to  failure.  This  has 
been  experimentally  identified  as  a  mode  of  failure  (reference  (h))  and  some 
investigators  (reference  (i))  have  applied  fracture  mechanics  methods  to 
quantitatively  define  failure  loads  for  various  angle-ply  laminates. 


COKCLUSIOKS  AND  RECOMMENDATIONS 

1.  finite  element  methods  can  be  used  in  ply-by-ply  analyses  of  local  edge 
effects  in  composite  laminates  and  x-rill  yield  valid  results.  In  the 
analysis  described  above,  only  the  simplest  one-  and  two-dimensional  con¬ 
stant-strain  elements  were  used,  Eve-n  better  results  should  be  obtainable 
with  siirpler  models  by  employing  newly-developed  parametric  solid  elements. 

2.  By  applying  supplementary  analytical  or  numerical  methods,  it  is  possible 
to  limit  tile  region  which  must  be  modelled  to  only  a  small  part  of  a  strtic- 
tural  element,  as  demonstrated  above.  Therefore,  a  finite-element  technique 
similar  to  that  used  in  the  analysis  of  the  straight  edge  is  potentially 
«ipplic^bi6  to  boundaries  of  arbitrary  geoiuetry* 

3.  The  edge  effect  stresses  resulting  from  the  analysis  presented  may  be 
sufficient  to  cause  delamination  and  subsequent  redistribution  of  stresses 
within  the  plies. 
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TABLE  I 

MATERIAL  PROPERTIES 


COMPOSITE 

FIBROUS  LAYER 

EPOXY  LAYER 

last 

20.00 

26.30 

0.607 

j  lElis  jL 

2.10 

2.56 

0.607 

tnsi 

2.10 

3.89 

0.607 

Gj^2* 

0.85 

1.05 

0.226 

Gj^3,  msi 

0.85 

1.05 

0.226 

G23,  nisi 

0.85 

1.05 

0.226 

V12 

0.21 

0.20 

0.340 

^13 

0.21 

0.15 

0.340 

Va3 

0.21 

0.11 

0.340 
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tabij:  II 

CONSTRAIKTS  APPLIED  TO  EDGE  EFFECT  MODEL 

CONDITION*  REGION  WHERE  APPLIED#- 


u=o 

=,  A  -  23h  ,,  -  - 

X  «  0;  yr  6h;  all 

o 

ii 

yz6h;  all  x,  z 

AM  -  0 

z  *  0;  all  X,  y 

Aii  p 

all  X,  y ,  z 

0  6  y  Ss  6h;  all  x,  z 

AX  ^C) 

AX  ^ 

X  =  Oj  y  »  all  z 

0  <  z  2h;  all  x,  y 

Ay 

X  »  0;  —  A  y  ^  6h;  1 

*u,  V,  w  are  displacements  in  x,  y,  z  coordinates  respectively 
+h  is  ply  thickness 
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FIGURE  4.  ELEMENT  ASSEMBLAGE  FOR  A  SINGLE  PLY 
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FIGURE  5.  ORTHOTROPIG  MATERIAL  STRESS-STRAIN  RELATIONS  (IN  MATERIAL  COORDINATES) 
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